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1 Introduction 



Electron-positron Bhahha scattering [Ij plays a special role in particle phenomenology. It 
provides a very efficient tool for luminosity determination at electron-positron colliders and 
thus it is crucial for extracting physics from the experimental data. Small angle Bhabha 
scattering is particularly effective as a luminosity monitor at the high energy colliders such 
as LEP and the future International Linear Collider (ILC) [2,3J. The large angle Bhabha 
scattering is used to measure the luminosity at colliders operating at a center of mass energy 

of a few GeV, such as BABAR/PEP-II, BELLE/KEKB, BES/BEPC, KLOE/DA<l>NE, 
and VEPP-2M [4J. The large angle Bhabha scattering will also be employed in order to 
disentangle the luminosity spectrum at ILC [5|6] . Since the accuracy of the theoretical eval- 
uation of the Bhabha cross section directly affects the luminosity determination, remarkable 
efforts have been devoted to the study of the radiative corrections to this process (see [2J for 
an extensive list of references) . QED contributions dominate the radiative corrections to the 
large angle scattering at intermediate energies 1-10 GeV and to the small angle scattering 
also at higher energies. The calculation of the QED radiative corrections to the Bhabha cross 
section is among the classical problems of perturbative quantum field theory and has a long 
history. The ffist-order corrections are well known (see [7j and references therein) . However, 
to match the impressive experimental accuracy reached at colliders, the complete second- 
order QED effects have to be taken into account. The evaluation of the two-loop virtual 
corrections is the main problem of the second-order analysis. The two-loop QED corrections 
can be divided into three main categories: (i) the corrections involving the electron vacuum 
polarization i.e. with at least one closed electron loop, (ii) the pure photonic corrections, and 
(iii) the corrections involving the vacuum polarization by heavy flavors (leptons or quarks) 
of mass rrif ^ rrie. Type (i) corrections have been evaluated in [8] including the full depen- 
dence on the electron mass mg. For the virtual corrections of type (ii) the full dependence 
on the electron mass is known with the exception of the two-loop box diagrams [11] Q At 
the same time, the complete result for the pure photonic corrections has been obtained in 
[T2] in the leading order of the small electron mass expansion through the infrared matching 
to the massless approximation [13] Jl It is sufficient for all the phenomenological applica- 
tions and has been recently conffimed within a slightly different framework [16]. For the 
heavy- flavor contribution the result is available only in the limit m'j -C s,t,u [TGlTTj . This 
condition, however, does not hold for r-lepton, c- and fe-quark in the practically interesting 
energy range of about a few GeV as well as for the top quark at the typical ILC energies 
500 GeV ^ ^ 1000 GeV. 

In this letter we consider the two-loop heavy-flavor contribution retaining the full de- 
pendence on m'j/s. In the next section we outline the structure of the corrections and the 
calculation method. In Sect. [3] we present the complete numerical result for the correction in 
two phenomeno logically interesting cases. Sect. H] contains our conclusions. The expansion 

''Partial results for the massive Bhabha scattering two- loop box diagrams are discussed in j9|10j . 

^The terms enhanced by logarithms of the electron mass have been derived in this approximation in 

mm- 
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of the result in the hmits mj <^ s and mj ^ s is given in the Appendix in analytic form. 

2 Structure of the Radiative Corrections and Calcula- 
tion Method 

We consider the phenomenologically interesting kinematical region s, t, m ~ mj ^ ml, 
where all the terms suppressed by the electron mass can be neglected. The perturbative 
expansion for the Bhabha cross section in the fine structure constant a is defined as follows 

where S^"'^ is the correction factor, with S^'^^ = 1. The leading order differential cross section 
in the small electron mass approximation takes the form 

where x = (1 — cos 9)/ 2 and 9 is the scattering angle. The virtual corrections are infrared 
divergent. These soft divergencies are canceled in the inclusive cross section when one adds 
the photonic bremsstrahlung [18]. The standard approach to deal with the bremsstrahlung 
is to split it into a soft part, which accounts for the emission of the photons with the energy 
below a given cutoff Scut ^ ^e, and a hard part corresponding to the emission of the 
photons with the energy larger than Scut- The infrared finite hard part is then computed 
numerically using Monte-Carlo methods with physical cuts dictated by the experimental 
setup. At the same time, the soft part is computed analytically and combined with the 
virtual corrections ensuring the cancellation of the soft divergencies in Eq. ([1]). We regulate 
all the soft divergencies by dimensional regularization m d = A — 2e space-time dimensions. 
The first-order heavy-flavor contribution to the cross section reads 

= ^ Ax{2x - l)Refn(/)(p)l + (2 - a;)Refn(;)(p/a;)l j , (3) 



where 



n?^(p) 



■- + Ap-{l + 2p) Jl-Ap \2 arctanhj ^ h i-n 

3 V V 1— 4p 




, (4) 



is the well known one-loop vacuum polarization function, p = mj/s, Qf is the heavy-flavor 
electric charge, and the number of colors is 1 for leptons and 3 for quarks. We work 
within the standard on-shell QED renormalization scheme. The second-order contribution 
can be split in the sum of two terms 

S^'' = € + 6i%^'^ , (5) 
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Figure 1: The two- loop diagrams associated with the logarithmic dependence of the corrections 
to the Bhabha scattering amplitude on me- Actually the diagram (a) is free of electron mass 
logarithms. The bold arrow circle corresponds to the heavy-flavor vacuum polarization. 

which correspond to the two-loop virtual correctioii^ and the one-loop correction to the single 
soft photon emission which factorizes into the product of the first-order contributions [T9] . 

The calculation of the virtual corrections is a highly nontrivial problem since it involves 
the two-loop four-point Feynman integrals that depend on four independent mass param- 
eters. However, the calculation is significantly simplified in the small electron mass limit. 
Eq. ([5]) can be rewritten according to the asymptotic dependence on rrie 



The logarithmic term in Eq. ([6]) is a remnant of the coUinear divergence which is regulated 
by the electron mass. The quantities and ^q^^ in Eq. ([6]) depend on s, t, and ruf only. 
The coUinear divergencies and hence the singular dependence of the corrections on me have 
a peculiar structure which was extensively studied in the context of QCD. In particular, the 
coUinear divergencies factorize and can be absorbed in the external field renormalization [20]. 
This property is crucial for our analysis because it allows us to perform the most difficult 
part of the calculation with a strictly massless electron. Indeed, due to the factorization, 
the singular dependence on nie is the same for the Bhabha amplitude and (the square of) 
the vector form factor [T^j. This attributes the total logarithmic corrections to the two- 
loop Bhabha scattering amplitude to the diagrams shown in Fig. [H However, due to the 
renormalization condition the vacuum polarization does not change the photon propagator 
near the mass shell where the coUinear divergencies are located. As a result, the diagram (a) 
is infrared finite even for rUe = and the singular terms are entirely contained in the one- 
particle reducible diagram (6). In calculating the cross section one must take into account also 
the contributions coming from the interference of the one-loop corrections to the amplitude 
and the corresponding soft emission; both contributions have a factorized form. Thus it is 
straightforward to obtain the coefficient of the logarithmic term in Eq. ([6]) which reads 






(7) 



'We do not consider the trivial correction given by two heavy-fermion loop insertions. 



3 



where e = At the same time the sum of the one-particle irreducible diagrams has a 

regular behavior in the small electron mass limit and can be computed with rrie = 0. Let 
us emphasize that this property in general holds only for the sum of the diagrams. The 
individual diagrams computed in a covariant gauge do exhibit the collinear divergencies for 
TTie = 0. This however does not pose any additional problem since we work in dimensional 
regularization. In this case the collinear divergencies show up as the extra poles in e which 
are not related to the soft emission and disappear in the sum of the one-particle irreducible 
diagrams. Thus all the "true" two-loop diagrams contribute only to the non logarithmic 
term in Eq. ([6]) which can be written as follows 

where the first term is determined by the soft emission and f{p,x) is a function of two 
dimensionless variables: p = m'j/s and x = —t/s. The two-loop problem with massless 
electron falls it the same complexity class as the one considered in Ref . [21] . The most difficult 
part of the calculation is represented, as expected, by the evaluation of the two-loop box 
graphs. By employing the Laporta Algorithm [22] to identify the set of master integrals (MI), 
it is possible to conclude that just two of the MI necessary for this calculation correspond 
to four point functions. These two are the only unknown ones among the Mis appearing 
in the calculation. They were evaluated by means of the differential equation method j23j 
and expressed in terms of Harmonic Polylogarithms [23] , and suitable generalizations of the 
latter [25]. The technical aspects of the calculation are going to be discussed in detail in 
|26j . Finally, the result for the function f{p,x) can be expressed in closed analytic form in 
terms of one- and two-dimensional generalized harmonic polylogarithms of maximum weight 
three. 

3 Results and Numerical Estimates 

The complete analytical result is rather lengthy and will be published elsewhere. In this letter 
we focus on the numerical impact of the correction in the phenomenologically important cases 
and we present the approximate analytical expression in different limits. 

As the first application, we consider the Bhabha scattering at the energy of 1 GeV which 
plays a crucial role in the determination of the hadronic vacuum polarization contribution to 
the muon anomalous magnetic moment [27p8] . The second-order contributions (a/7r)^5*-^'' 
of heavy leptons and quarks to the differential cross section are plotted as function of the 
scattering angle in Figs. [2] and [3] for ^/s = 1 GeV, which is the center-of-mass energy 
of the KLOE experiment at DA$NE. All the terms involving a power of the logarithm 
In^Ecut/^) are excluded from the numerical estimates because the corresponding contribution 
critically depends on the event selection algorithm and cannot be unambiguously estimated 
without imposing specific cuts on the photon bremsstrahlung. Note that the energy under 
consideration is sufficiently below the quarkonium threshold and the heavy-quarks can be 
treated perturbatively. The contributions of the r-lepton, c and b quarks are suppressed 
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Figure 2: Two-loop corrections to the Bhabha scattering differential cross section at y/s = 1 GeV 
due to a closed loop of r-lepton (dotted line), c-quark (dashed line) and 6-quark (solid line) for 
rUc = 1.25 GeV and rub = 4.7 GeV. 
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Figure 3: Two-loop corrections to the Bhabha scattering differential cross section at ^/s = 1 GeV 
due to a closed loop of muon (dashed line). The solid line represents the sum of the contributions 
of the muon, r-lepton, c-quark and 6-quark. 
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Figure 4: Two-loop leptonic corrections to the Bhabha scattering differential cross section at ^/s = 
500 GeV. The dash-dotted line represents the electron contribution including the soft-pair radiation. 
The dashed and dotted lines represent the contributions of muon and r-lepton, respectively. The 
solid line is the sum of the three contributions. 
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Figure 6: Two-loop corrections to the Bliabha scattering differential cross section at 6 = 60° due 
to a closed loop of muon. The solid line represents the exact result. The dashed and dotted lines 
represent the results of the large-mass expansion (Eq. (fT3]) ) and small- mass expansion (Eq. (fTOll ). 
respectively. 



with respect to the muon at least by one order of magnitude and they are given separately 
in Fig. [2j Thus the total heavy-flavor contribution is dominated by the muon loop and it 
reaches 0.45 permille in magnitude at 6* ~ 140°, see Fig. |3l 

Let us now discuss Bhabha scattering at high energies characteristic to the ILC. We 
consider only the contributions of the leptons and the top quark because the lighter quarks 
need a special treatment due to hadronization effects [29] . In Fig. H] we present the plots for 
the contributions of muon, r-lepton, and electron (including the soft electron-positron pair 
emission) at ^/s = 500 GeV. The total leptonic correction reaches 1.3% at 6* ~ 120° and it 
is comparable with the pure photonic term [12]. The contribution of top quark is plotted in 
Fig. [5l It is significantly smaller than the leptonic one. 

The analytical structure of the result becomes much simpler in the small-mass limit 
m'j <^ s and the large-mass limit mj ^ s. The corresponding expansions of the function 
f{p,x) at p —>■ and p —>■ oo are given in the Appendix. Both expansions break down 
near the threshold ^/s ~ 2m f. In Fig. [6] the muon contribution at the fixed scattering angle 
9 = 60° is plotted as function of in the threshold region around 2m^ ~ 216 MeV. In this 
plot we present the exact result along with the first nonvanishing terms of the small and large 
mass expansions, Eqs. (fTOl [T3i) . We observe that the asymptotic results fail to approximate 
the exact one with 10% accuracy in the interval 1.7 ^ y/s/nif <, 10. In particular they are 
completely useless for the description of the top quark contribution at the typical energies 
of the ILC. 
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4 Summary 



To conclude, we have derived the two-loop radiative corrections to Bhabha scattering due 
to the heavy-flavor vacuum polarization in the leading order of the small electron mass 
expansion and for any ratio of the heavy-fermion mass to the energy of the process. This 
completes the calculation of the QED part of the two-loop corrections. We have analyzed 
the phenomenological impact of our result in two different experimental setups. For e'^e~ 
colliders operating at the energies of a few GeV, the correction is dominated by the virtual 
muon loop and reaches 0.45 permille at large scattering angles. At the typical ILC energies 
the r-lepton and muon loops result in the comparable contributions of about two permille 
so that the total leptonic correction amounts up to 1.3%. The contribution of the top quark 
is smaller and at y/s = 500 GeV reaches approximately 0.25 permille in magnitude. 
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Appendix 

The small-mass expansion of the function f{p,x) is of the following form 

oo 
n=0 

where fn{p,x) depend on p only logarithmically. For the leading term we obtain 
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where Q/ is the charge of the heavy fermion in units of the electron charge. Eq. flTUl) is in 
agreement with the result of Refs. [TBIT7] . The next-to-leading term is new and reads 
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The expansion in the large-mass hmit takes the form 

oo 

f{p,x) = Y.p--Up,x), (12) 

n=0 

where the leading n = term vanishes because of the renormalization condition and x) 
depend on p only logarithmically. For the next-to-leading term we obtain 
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